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Abstract 

Most experimental data on the single-molecule magnet Mni2 acetate have been successfully 
explained by the assumption that the Mni2 acetate has an effective ground-state spin of S" = 10. 
However, the effect of the low-lying excited manifolds caused by interactions between Mn spins 
has not been well understood. To investigate the features of the low-lying excited manifolds, the 
intramolecular exchange interactions are calculated using density- functional theory (DFT). With 
the calculated exchange parameters, the energy gap between the = 10 ground-state and the 
first excited-state manifold is calculated by diagonalization of the Heisenberg Hamiltonian. The 
upper limit on the energy gap is about 40.5 K which is likely to be overestimated due to incomplete 
treatment of the Coulomb potential within DFT. It is found that there are several 5 = 9 low-energy 
excited-state manifolds above the 5 = 10 ground-state manifold. The magnetic anisotropy barriers 
for the low- lying spin excitations are calculated using DFT. Based on the calculations, it is found 
that the anisotropy barriers for the low-lying excited manifolds are approximately the same as that 
for the ground-state manifold, which is applicable to other single- molecule magnets such as Mn4. 

PACS numbers: 75.50.Xx, 71.15.Mb, 75.30.Gw, 75.30.Et 
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I. INTRODUCTION 



Single-molecule magnets (SMMs) have been extensively studied for the past decade 
because of both scientific and practical reasons: macroscopic quantum phenomenal^ 
and possible utilization as magnetic storage devices.- A prototype of the SMMs is 
[Mni20i2(CH3COO)i6(H20)4] ■2(CH3COOH)-4(H20) (hereafter Mni2),-^ which is a three- 
dimensional array of identical S* = 10 molecules. Very recently, derivatized Mni2 type 
molecules were successfully deposited on a gold film,^ which further enhances the prospects 
for storing magnetic information in a single molecule. A single molecule of Mni2 has four 
ferromagnetically coupled Mn^"*" ions (5 = 3/2) in the cubane and eight ferromagnetically 
coupled Mn^+ ions {S = 2) in the crown as schematically shown in Fig. ^ In the ground 
state, the magnetic moments of the eight Mn^+ ions are antiparallel to those of the four 
Mn*^+ ions which leads to a total spin of S* = 10. 

So far many interesting features of the SMM Mni2 such as magnetization steps in the 
hysteresis loops,"^ electron paramagnetic resonance (EPR) transitions,'''^ and low-energy 
excitations in inelastic neutron scattering ,^i^°i^^ have been well understood by considering 
each molecule as an object with an effective spin oi S = 10. However, magnetic susceptibility 
measurement o^^i^^i^"^ have demonstrated that some experimental data could not be in accord 
with the 5* = 10 single-spin picture and suggested that the first excited-state manifold may 
be located at 35 K— 40 K above the S" = 10 ground-state manifold within which the energy 
barrier to magnetization reversal is 65 K in zero field. Since the first excited-state manifold 
overlaps with the ground-state manifold above 35 K— 40 K, an internal many-spin structure 
of a single molecule should be included to explain the high-energy experimental data. The 
dimension of the Hilbert space of this SMM is so large that a simplified eight-spin model 
with strong Dzyaloshinsky-Morya interaction was first proposed to include many-spin effects 
by dimerizing two strongly bonded Mn spins.J- Although this eight-spin model explained 
some experimental data, it was limited to features below 50 K. At higher temperatures the 
dimerization scheme breaks down. Later all twelve Mn spins were included in the Heisenberg 
Hamiltonian and excited-state manifolds were clarified by diagonalization of the Hamiltonian 
with constraining the fixed energy gap between the first excited-state and the ground-state 
manifold of 35 Nonetheless, there is still a big controversy over the energy gap between 

the first excited-state and the ground-state manifold, and there has been less exploration 
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of high-energy features caused by many-spin effects. Here "high-energy" means low-lying 
excited manifolds above the = 10 ground-state manifold. Recent high-energy inelastic 
neutron scattering measurements exhibited a broad anomalous peak at 10 K, which may 
have a magnetic origin (not due to phonons) and could not be rationalized by the S* = 10 
single-spin picture because the first allowed peak by the picture is about 14 K.^ More recent 
high-field EPR measurements revealed that several transitions could not be justified by the 
S* = 10 single-spin manifold but that they were well explained if we assumed the S = 9 
first excited-state manifold to be situated at 10-16 K above the ground-state manifold.— 
Additionally, there have been many attempts to determine the signs and magnitudes of the 
exchange couplings between Mn ions from different approachesi^ii^iiSiSS but no consensus 
has arisen. 

To elucidate the energy gap and examine properties within the low-energy excited mani- 
folds, a single molecule of the Mni2 is considered and the exchange interactions between Mn 
spins are calculated using density-functional theory (DFT), which are presented in Sees. II 
and III. Our calculated values are compared with those from other groups. Through 
diagonalization of the isotropic Heisenberg exchange Hamiltonian with our calculated ex- 
change constants using the Lanczos method , ^^^^^ excited-state manifolds are manifested and 
energy gaps between different manifolds are calculated in Sec. IV. The magnetic anisotropy 
barriers for low-energy spin excitations are calculated using DFT and the anisotropy barri- 
ers for low-lying excited manifolds are discussed in Sec. V. Our conclusions are presented in 
Sec. VI. 

II. DFT CALCULATIONS 

Our DFT calculations^^, are performed with spin-polarized all-electron Gaussian-orbital- 
based Naval Research Laboratory Molecular Orbital Library (NRLMOL)24 which is ideal for 
studying a single molecule to a small number of unit cells. Here we use the Perdew-Burke- 
Ernzerhof (PBE) generalized-gradient approximation (GGA) in the exchange-correlation 
potential.— Since the SMM Mni2 has fourfold symmetry, for fully symmetrized calcula- 
tions the total number of inequivalent atoms to consider is reduced to 176/4 = 44. To 
save geometry-optimization time without losing interesting physical properties, the follow- 
ing simplified form of the SMM Mni2 is used: [Mni20i2(HCOO)i6(H20)4] (acetic acids and 
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water molecules of crystallization are not included, and 16 acetates, CH3COO, are replaced 
by 16 formates, HCOO). The zero-field total anisotropy barrier for the = 10 ground-state 
manifold does not change much with this simplification. The total magnetic moment for the 
ground state was confirmed to be 20yUB, which is in good agreement with experiment and 
it is stable. Details of the optimization schemes and electronic properties of the optimized 
geometry for the ground state were discussed elsewhere.— Hereafter, unless specified, our 
calculations have been carried out for the above simplified form. 

III. INTRAMOLECULAR EXCHANGE INTERACTIONS 

The fourfold symmetry and the geometry of a single Mni2 molecule indicate that there are 
three symmetrically inequivalent Mn sites and four different exchange interactions between 
Mn spins as shown in Fig. ^ To calculate the exchange constants using DFT, it is assumed 
that the magnetic moments of all Mn ions are aligned along a particular direction (collinear). 
Then eleven distinctive spin configurations are constructed by reversing magnetic moments 
of a few Mn ions simultaneously from the 5 = 10 ground-state. (See Table |l} For example, a 
Ms = 9 state {Mg is an eigenvalue of the total spin operator 5* projected along the particular 
direction) can be built by flipping both one Mn spin in the cubane and another Mn spin in 
the crown within the S = 10 ground state: = -3/2 x (3 - 1) + 2 x (7 - 1) = 9. There 
are several ways to construct Mg = 9 states with collinear spins that are labeled as 9-b, 9-c, 
and 9-d in Table HI All of the examined collinear spin configurations are not eigenstates of 
5*^. Mean-field calculations are not directly applicable to optimizing excited states such as 
15 = 9, Ms = 9) in Mni2, because those eigenstates are represented by linear combinations of 
many single Slater determinants. The geometries for the distinctive spin configurations are 
taken to be the same as that for the ground state except for the proper spin arrangements. 
It is confirmed that a slightly different initial geometry for the collinear spin configurations 
does not siginificantly alter our calculated values of the exchange constants. Since fourfold 
symmetry is broken for the spin configurations considered, unsymmetrized calculations are 
required with a total number of 100 atoms. The energies of the spin configurations are self- 
consistently minimized with a small basis set and fine mesh.™ Our calculations show that 
the minimized energies range from 0.04 eV to 0.2 eV above the ground-state ferrimagnetic 
structure. From Ref. it is known that for the 5 = 10 ground state the lowest electronic 
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excitation is a majority spin excitation between the eg levels in Mn(3d) states, which is 
about 0.44 eV (majority LUMO-HOMO gap), and that the energy gap between minority 
lowest unoccupied molecular orbital (LUMO) and majority highest occupied molecular or- 
bital (HOMO) is ~ 0.89 eV. Comparison of the minimized energies to the LUMO-HOMO 
gaps indicates that these different spin configurations are lower-energy spin excitations than 
moving an electron from the majority HOMO to the unoccupied majority (or minority) 
orbital. Thus they will be called low-energy spin excitations. There are twelve equations to 
solve for five unknowns (the background energy, Eq, four exchange constants Ji, J2, J3, and 
J4) so a least-square- fit (LSF) method^^ is used. 

Our calculated values of J's are shown in Table [H] in comparison with results of other 
groups. When the minimized energies of the low-energy spin excitations are recalculated 
using the calculated values of J's, it is found that they are in good agreement with the 
DFT-calculated energies as indicated in Table IH Prior to comparison to the results of other 
groups, let us briefly review the essence of the different approaches from the different groups. 
In Ref. Q the value of Ji was first determined from the experimental data, and then Jo and 
J3 were varied with fixed value of J4 = to reproduce the S* = 10 ground state. In Ref.llfilthe 
values of J's were slightly varied from the previously reported values^^ using diagonalization 
of Heisenberg exchange Hamiltonian to provide the = 10 ground state and the energy gap 
of 35 K between the ground-state and the S = 9 first excited-state manifold. In Ref. |l7| the 
exchange parameters were obtained to reproduce their megagauss experimental data with a 
constraint of the energy gap of 35 K. In Ref. the LDA+U method was used to include 
electron correlations between Mn ions, and the values of the intra-atomic Hund's exchange 
parameter J and the average Coulomb parameter U were determined to obtain the correct 
magnetic moment for the 5* = 10 ground state. Then the exchange constants were calculated 
from the values of J and U. Our calculations show that the coupling between one Mn spin 
in the cubane and another Mn spin in the crown closest to the Mn spin in the cubane, Ji, is 
confirmed to be the strongest antiferromagnetic coupling, although its absolute magnitude 
does not agree among different groups. Our calculated value of Ji agrees somewhat with 
those from Refs. Il7ll20l . The coupling between a Mn spin in the cubane and the other type 
of Mn spin in the crown, J2, is confirmed to be weaker antiferromagnetic than Ji. But 
the rate of J2 to Ji does not agree among different groups. Our value of J2 is between 



the values obtained from Refs. 
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121120. The coupling between Mn spins in the cubane, J3, 
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turns out to be the weakest and is ferromagnetic. The couphng between Mn spins in the 
crown, J4, is weakly antiferromagnetic. Overall, the DFT calculated values are closest to the 
results obtained from Ref. 0, and they are 20%-40% larger than those from the LDA+U 
method except for the weakest Ja- However, they are significantly smaller than the other 
two reported results from Refs. I16II191 Our calculated values of J' may be overestimated 
compared to experimental data for the following reasons. It has been found from other SMMs 
such as Mu/i [Ref. and Vi5[Ref. |3fl] that the DFT-calculated exchange interactions were 
overestimated by a factor of 2 and 3 compared to the experimental data respectively.—!^ 
This overestimated exchange interactions are due to the fact that the PBE generalized- 
gradient approximation does not fully treat self-interaction corrections of the localized d 
states in the Coulomb potential. Proper treatment of the electron correlations in Mn atoms 
may improve our DFT calculations of the exchange interactions. This is supported by the 
fact that the LDA+U method produces 20%-40% smaller exchange interactions than our 
DFT calculations and that the exchange interactions decrease with increasing the value of 
U in the LDA+U method.^ 



IV. HEISENBERG EXCHANGE HAMILTONIAN 



To calculate the energy gap between the ground-state and the first excited manifold, we 
use the following isotropic (i.e. no anisotropy) Heisenberg exchange Hamiltonian: 

^ = JijSi ■ (1) 

i<j 

where the sum runs over all pairs connected by the exchange interactions and the values of 
the coupling constants Jij labeled in Fig^ are in Table IHl The exchange Hamiltonian is a 
10® X 10® matrix, which is too large to diagonalize using currently available computers. Thus, 
using the fact that is a good quantum number, we classify the total number of the 10® 
states into 45 different constant Mg states such as Mg = ±22, ±21, ±20, ±10, ±9, ±1, 0. 
Since we are interested in low-lying excited manifolds, we examine Mg = —12, Mg = —11, 
Mg = —10, Mg = —9, and Mg = —8 states only. If the 5* = 10 manifold is the ground state, 
then the lowest-energy eigenstate of the exchange Hamiltonian with a basis set of all Mg = 
—10 states corresponds to \S = 10, = —10). Additionally the lowest-energy eigenstate 
of the Hamiltonian with a basis set of all Mg = —9 states corresponds to 15* = 10, Mg = —9) 
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and it should have the same energy as IS* = 10, Ms = —10). The eigenstates with the same 
S must be degenerate if they belong to the same manifold. The first-excited eigenstate of 
the Hamiltonian with the basis set of all Mg = —9 states corresponds to IS* = 9, = —9) 
if the S = 9 manifold is the first-excited manifold above the S* = 10 ground state. There 
are a total number of 269148 = 12 states, 484144 = -11 states, 817176 = -10 
states, and 1299632 M, = -9 states, 1954108 M, = -8 states. 

We diagonalize the full Hamiltonian matrices for the Mg = —12, Mg = —11, Mg = —10, 
Mg = —9, and Mg = —8 states using the Lanczos method . A few calculated low-lying 
energy eigenvalues are shown in Table ITTTl The lowest energy for the Mg = — 10 states is 
identical to those for the Mg = —9 and Mg = —8 states, but it is much lower than the lowest 
energies for the Mg = — 11 and Mg = —12 states. The lowest energies for the Mg < —12 
states will increase with increasing \Ms\. This indicates that the ground-state has S = 10. 
Next higher energies are from the Mg = —9 and Mg = —8 states. The three excited energy 
eigenvalues for the Mg = —9 states coincide with those for the Mg = —8 states which are 
much lower than the excited energies for the Mg = —10, Mg = —11, and Mg = —12 states. 
This suggests that the three excited states for the Mg = —9 {Mg = —8) states belong to 
S = 9 manifolds. So the first excited-state manifold has S = 9 located at 40.5 K above the 
S" = 10 ground state, and there are two more S = 9 excited manifolds located at 7 K and 
40 K above the S = 9 first excited manifold (Fig. Our calculated values of the energy 
gaps are upper limits because DFT calculations may overestimate the exchange interactions 
due to the reasons discussed in Sec. III. If the exchange constants are reduced by a half, 
then the energy gaps will be also reduced by a half. Single-ion anisotropy parameters which 
we did not include in Eq. does not significantly affect the energy gap. 

V. MAGNETIC ANISOTROPY BARRIER 

The magnetic anisotropy barriers for the eleven low-energy spin excitations as well as the 
ground state are calculated with the assumption that the barriers are caused by the spin- 
orbit interaction only. Other effects such as noncollinearity of the magnetic moments of Mn 
ions and spin-orbit-vibron coupling^^ on the anisotropy barrier will not be included in our 
study. However correlation effects due to the addition of multideterminantal wavefunctions 
are addressed in part by the further ongoing discussion. The spin-orbit interaction Vls can 
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be decomposed into one-electron operators fi related to the electric fields caused by nuclei 
and two-electron operators gij related to the electric fields due to the rest of electrons as 
follows: 

Vls = -^T.S^■{p^y<^<^^) (2) 

= T.f^ + T.9^j, (3) 

^•^ = E^^ + Er3^ (4) 

U = -E^^^-fv^x ^'-'^A (5) 

where the summation for i (z/) runs over all electrons located at (all nuclei located at Ry 
with nuclear number Z^,), c is the speed of light, Si is the spin operator of the iih electron, and 
Pi is the momentum of the zth electron. Within the self-consistent field (SCF) approximation, 
we have calculated the wave function ^i^, = ipi„(r)xa satisfying TCl'^ia) = ejo-I^Pjo-), where H 
is the many-electron Hamiltonian without spin-orbit coupling, ipi^ is a spatial function and 
Xa is a spinor. Therefore, the energy shift, A2, due to Vls using second-order perturbation 
theory is written by2& 

_ ^ ^ {^^. I Vls I "^ja' ) i^ja' \ Vls I ^i.) 

where ejo- is the energy of the occupied state with spinor a, and ej„' is the energy of the 
unoccupied state with spinor a'. The summation runs over all occupied and unoccupied 
states of all atoms within a certain energy window and over up and down spinor states. The 
contributions for the cases that both of i and j are occupied (or unoccupied) to the energy 
shift cancel each other out. For uniaxial systems such as Mni2 the second-order energy shift 
is simplified as follows: 

A2 = D{S,f (8) 

where D is the uniaxial anisotropy parameter which can be calculated from Eq. (0). For the 
= 10 ground state, as shown in Table H] the uniaxial parameter is about 0.54 K, which is 
in good agreement with experiment.— 
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An eigenstate |$) of the many-electron Hamiltonian that belongs to an excited-state 
manifold can be written in terms of a linear combination of many Slater determinants: 
1$) = J^^C^l^^), where |$^) is a single Slater determinant (orthonormal basis function), 
and some of |$^) correspond to our examined low-energy spin excitations. Therefore, the 
expectation value of the spin-orbit coupling Vis with respect to 1$), is decomposed into 
diagonal and off-diagonal elements such as 

{<^>\Vls\<i>) = E c;c,{<^,\Vls\'^,) + E Vl5|$.) , (9) 

where the summations run over spatial coordinates and spin variables of all electrons in 
Mni2. Vls is made of one-electron and two-electron operators as shown in Eq. 0. To 
facilitate calculations of the diagonal and off-diagonal elements, the off-diagonal elements 
of the one-electron operator, J2i{^fi\fi\^u)^ are explicitly written in terms of single-electron 
spin-orbit als.— 



''"ir--jO'JV 




(10) 



where Vj is also applied to the single-electron spin-orbitals and the normalization factor 
1/A^! cancels out A^! permutations of the single-electron spin-orbitals in ($^|. As can be 
seen from Eq. (HU)), because all spin-orbitals V'laiXai,- ■ ■ V^TVa^vXajv; ^U^^i''"' V^W^U' 
are orthonormal to each other, the off-diagonal elements of the one-electron operator /j 
are nonvanishing when only one single-electron spin-orbital is different between |$^) and 
\^v)- The off-diagonal elements of the two-electron operator gij are nonvanishing when no 
more than two single-electron spin-orbitals are different between |$^) and Different 
single Slater determinants (or different low-energy spin excitations) are not connected by 
the one- or two-electron operators, so {^^\Vls\^v) = for /i 7^ z/. Hence, the off-diagonal 
contributions to the energy shift A2 is zero. As shown in Table |lj the calculated magnetic 
anisotropy barriers for the low-energy spin excitations are almost the same as that for the 
ground state 5* = 10 manifold. This indicates that the diagonal elements of the spin- 
orbit interaction [Eq. 0] do not vary much with different spin excitations: ($^|Vl5|$^) ~ 
{^u\Vls\^,) for Then {^Vls\^) ~ ($^| V^^l^^)© JC'^H where EJC^mI' = 1 by 
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definition. So the magnetic anisotropy barriers for low-lying excited-state manifolds must be 
approximately the same as that for the ground-state manifold. Note that we have calculated 
second-order anisotropy barriers only so that the barriers shown in Table E] are all between 
54 K and 55 K. The results of this discussion suggest that the DFT-based determination of 
MAE is valid under the following conditions: (i) single-electron and two-electron excitations 
cost an order of eV, (ii) considered collinear spin excitations cost an order of magnitude less 
energy than the single-electron and two-electron excitations. 



VI. CONCLUSION 



We have calculated the intramolecular exchange interactions among Mn spins in the 
single-molecule magnet Mni2 ([Mni20i2(HCOO)i6(H20)4]) considering localized low-energy 



spin excitations using density-functional t 
constants agree with the results from Refs. 



leory. Our calculated values of the exchange 



2(1 . With our calculated exchange constants, 



we have diagonalized the isotropic Heisenberg Hamiltonian using the Lanczos method. We 
have confirmed that the ground state is S = 10 and that there is a S* = 9 first-excited 
manifold located at 40.5 K above the 5 = 10 ground-state manifold. This energy gap 
(40.5 K) is an upper limit because it is likely reduced by including more exact treatments of 
the Coulomb potential or electron correlations. This indicates that the first excited manifold 
may be situated much lower than the high energy levels within the 5 = 10 ground-state 
manifold. We have also calculated the magnetic anisotropy barriers for the low-energy 
spin excitations using DFT and the second-order perturbation theory with the assumption 
that the spin-orbit coupling is the most important interaction. Our DFT calculations have 
showed that the anisotropy barriers for the low-energy spin excitations are approximately 
the same as that for the S" = 10 ground state. From this result and our understanding of 
the spin-orbit interaction, we conclude that the anisotropy barriers for the low-lying excited 
manifolds are approximately same as that for the S* = 10 ground state. 
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TABLE I: Here Mg = 6 — c denotes Mj, = 6, where Mg is an eigenvalue of the total spin operator 
projected along the easy axis and "c" is attached to distinguish between different types of Mg = 6 
states. From the second column, shown are flipped spins to create low-energy spin excitations 
relative to the S" = 10 ground-state labeled in Fig^ Ising energy expressions, DFT-calculated 
energies relative to the ground state, AE (in units of eV), the energy differences between DFT 
results and least-square-fit (LSF), E^^"^— E^"^^, and the magnetic anisotropy barriers (MAE) (in 
units of K). 
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-0.0049 


54.1 


6-c 


9 


Eo 


- 12 Ji - 48 J2 + 27J3 + 32 J4 


0.0352 


-0.0064 


54.8 


9-b 


1, 9 




Eq - 24Ji - 24J2 + 36 J4 


0.0598 


-0.0134 


54.5 


6-b 


5 


Eo 


- 24 Ji - 24J2 + 27J3 + 32 J4 


0.0780 


0.0095 


55.2 


5-b 


1, 5, 9 




Eo - 24Ji - 24J2 + 32 J4 


0.0902 


0.0171 


54.8 


5- a 


1, 7, 9 




Eo - 24 Ji 


0.1328 


-0.0045 


55.4 


8-b 


1, 4, 5, 9 


Eq 


- 12 Ji - 24 J2 - 9 Js + 36J4 


0.1361 


0.0014 


54.1 


8 


1, 4, 9, 12 




Eo - 24Ji - 9 J3 


0.1377 


-0.0012 


54.9 


9-c 


4, 5 




Eo - 12 Ji - 24J2 + 36J4 


0.1445 


0.0114 


54.9 


13 


1 




E0-I2J1- 24 J2 + 64 J4 


0.1496 


-0.0064 


53.6 


8-c 


1, 8, 9 




Eo - UJi - 9J3 


0.1929 


-0.0057 


54.8 


9-d 


4, 9 




Eo - 24J2 + 32 J4 


0.1955 


0.0027 


55.0 
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TABLE II: Comparison of our DFT-calculated intramolecular exchange couplings with those from 
LDA+U calculations^'^, diagonalization of the Heisenberg mode^^^, and experimental resultsi^ 
in units of K. A positive sign denotes an antiferromagnetic coupling. 





DFT 


LDA+U[Ref.2^ 


Ref.JJ 


Ref.J^ 


Exp.[Ref.i9| 


Jl 


115 


94 


119 


430 


432 


J2 


84 


52 


118 


170 


173 


•h 


-4 


60 


-8 


170 


166 



J A 17 14 23 -129 



TABLE III: Several lowest-energy eigenvalues of the Hamiltonians for a total number of N Mg = 



-12, M, = 


-11, Ms = -10, Ms = -9, and Ms = 


—8 states. 


Ms 


N 


Energy (K) 


-12 


269148 


-2574.2988 


-11 


484144 


-3034.6061 






-2926.5083 


-10 


817176 


-3464.9569 






-3034.6061 


-9 


1299632 


-3464.9569 






-3424.4283 






-3417.4348 






-3384.2171 


-8 


1954108 


-3464.9569 



-3424.4283 
-3417.4348 
-3384.2171 
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FIG. 1: Schematic diagram of the exchange interactions within a single molecule for the SMM 
Mni2. The empty circles denote Mn^+ ions, the filled circles denote Mn^+ ions, and each Mn ion 
is numerically labeled. The bond length between spin 1 and 9 is shorter than that between spin 1 
and 6. The thick solid lines are for Ji, the thick dashed lines are for J2, the thin solid lines in the 
inner cubane are for J3, and the thin solid lines in the outer crown are for J4. 





33K 




r 





S=9 



S=9 
S=9 



40.5K 



* S=10 

FIG. 2: Calculated low-lying excited manifolds for the SMM Mni2. 
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